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We report on theoretical Auger electron kinetic energy distribution originated from 
sequential two-step Auger decays of molecular double core-hole (DCH) state, us- 
ing CH4, NH3 and H2CO molecules as representative examples. For CH4 and NH3 
molecules, the DCH state has an empty Is inner-shell orbital and its Auger spectrum 
has two well separated components. One is originated from the 1st Auger transition 
from the DCH state to the triply ionized states with one core hole and two valence 
holes (CVV states) and the other is originated from the 2nd Auger transition from 
the CVV states to quadruply valence ionized (WW) states. Our result on the NH3 
Auger spectrum is consistent with the experimental spectrum of the DCH Auger 
decay observed recently [Phys. Rev. Lett. 105, 213005 (2010)]. In contrast to CH4 
and NH3 molecules, H2CO has four different DCH states with Cls"^, 01s~^ and 
Cls"^01s~^ (singlet and triplet) configurations, and its Auger spectrum has more 
complicated structure compared to the Auger spectra of CH4 and NH3 molecules. 
In the H2CO Auger spectra, the Cls^^Ols^^ DCH — t- CVV Auger spectrum and 
the CVV — 7- WW Auger spectrum overlap each other, which suggests that iso- 
lation of these Auger components may be difficult in experiment. The Cls^^ and 
Ols^^ DCH —7- CVV Auger components are separated from the other components 
in the H2CO Auger spectra, and can be observed in experiment. Two-dimensional 
Auger spectrum, representing a probability of finding two Auger electrons at specific 
pair of energies, may be obtained by four-electron coincidence detection technique 
in experiment. Our calculation shows that this two-dimensional spectrum is useful 
in understanding contributions of CVV and WW states to the Auger decay of 
molecular DCH states. 
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I. INTRODUCTION 



Molecular double core-hole (DCH) state is a state of a molecule with two inner-shell 
vacancies; they mostly correspond to two K-shell vacancies in the literature,-"— although 
two L-shell vacancies or one K-shell and one L-shell vacancies have been treated as well.— 
Two different kinds of molecular DCH state exist: single-site (ss) DCH state having core-hole 
vacancies at the same atomic site and two-site (ts) DCH state with two core-hole vacancies 
at different atomic sites. Properties of molecular DCH state with K-shell vacancies was first 
studied by Cederbaum et al.J. at 1986 and have been discussed occasionally since then.-"- In 
contrast, experimental realization of probing molecular DCH state has not been possible until 
recently because of technological difficulties. Recently, X-ray free electron laser (XFEL) at 
Linac Coherent Light Source (LCLS) started operatioUj^- which has capability to generate 
high intensity, short laser pulses required to produce molecular DCH state by sequential 
two-photon two-electron ionization. Using this XFEL facility. Fang et al.- identified the 
ss-DCH state of nitrogen molecule by analyzing photoelectron and Auger electron spectra. 
Cryan et al.^ also obtained the Auger electron angular distribution originated from the 
N2 DCH Auger decay in the molecular frame. At the same time, it has become possible 
to study molecular DCH state using synchrotron radiation (SR), combined with multi- 
electron coincidence technique. In contrast to two-photon sequential ionization by XFEL, 
SR mainly creates molecular ss-DCH state through single-photon two-electron ionization 
process. Eland et al.— detected the ss-DCH states of NH3 and CH4 molecules in their 
SR experiment. For NH3 molecule, they obtained the Auger electron spectrum originated 
from the cascade Auger decay of the ss-DCH state. The 2h-lp pre-edge resonance state of 
NH3, with double core- hole single electron valence excited configuration, was also identified, 
where they found that the state decays predominantly through the spectator Auger process. 
In different SR experiment, Lablanquie et al.— determined binding energies of the ss-DCH 
states for N2, O2, CO and CO2 molecules. Using the four-electron coincidence detection 
method, with two photoelectrons and two Auger electrons, they obtained two-dimensional 
(2D) Auger intensity distribution of the N2 ss-DCH decay as functions of two Auger electron 
kinetic energies. Their photoelectron spectrum of N2 shows clear signature of the DCH 
shake-up satellite states. Based on the Auger spectrum of the N2 DCH satellite state, they 
concluded that the satellite state decays by the spectator Auger process. Although number 
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of experimental studies on molecular DCH state is still limited, it is expected to increase 
in near future because XFEL facilities are being constructed in Japan and EU/^ and also 
because more SR experiment will be performed to study molecular DCH state. 

Several theoretical studiesi^-ii^ have been performed on molecular DCH state to clarify its 
basic properties. As in single core-hole state, molecular DCH state is not stable but decays 
by Auger electron emission or X-ray emission, where Auger decay is expected to be dominant 
for low-Z elements.— At present, number of theoretical works on Auger decay of molecular 
DCH state is quite limited, where only two studies exist: theoretical N2 DCH Auger spectrum 
by Fang et al.* and characterization of transient and final ion states in CH4 and NH3 DCH 
Auger decays by Eland et al.— Detailed knowledge about Auger electron kinetic energy 
distribution, such as assignment of peaks and prediction of Auger intensity, is important 
to interpret experimental molecular DCH Auger spectrum, in which single core-hole Auger 
spectrum or other instrumental noise often prevents straitforward extraction of DCH Auger 
spectrum.— In addition, detailed theoretical information may be valuable for studying time- 
dependent process of molecular DCH formation and decay in XFEL experiment. 

Auger decay of molecular DCH state mainly involves three different kinds of electronic 
states:—"^ DCH state, triply ionized states with one core hole and two valence holes (CVV 
states) and quadruply valence ionized (WW) states. These electronic states and their 
relations are schematically shown in Fig. [1] In the first Auger transition, the DCH state 
decays to the CVV states, then these CVV states decay to the WW states in the second 
Auger transition. Thus, two Auger electrons are emitted in molecular DCH Auger decay, 
in contrast to one Auger electron emission in Auger decay of single core-hole state. Other 
processes may also be involved in molecular DCH decay, such as X-ray emission, double 
Auger process^'' where a valence electron fills the vacant core-hole with emission of two 
Auger electrons, or direct Auger decay from DCH state to WW state where two valence 
electrons simultaneously fill the vacant core-holes with an emission of an Auger electron. 
These processes are not treated in the present work, because the probability of X-ray emission 
is expected to be small for low-Z elements treated in the present work. Aso, description of 
the double Auger process and the direct Auger decay from DCH state to WW state 
requires the second order perturbation in terms of the Coulomb interaction,— and thus their 
transition probabilities are expected to be small compared to normal Auger process which 
can be described by the first order perturbation. 
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In this work, Auger electron kinetic energy distributions are calculated and analyzed for 
the DCH states with K-shell vacancies of CH4, NH3 and H2CO molecules. Auger decay 
in CH4 and NH3 involves only ss-DCH state. Auger decay in H2CO, on the other hand, 
involves both ss- and ts-DCH states and therefore more complex Auger spectrum is expected. 
Although number of molecules is limited, our results may provide useful insights to discuss 
Auger decays of DCH states in other molecules. 

We employ the complete active space self-consistent field (CASSCF) and configura- 
tion interaction (CASCI) methods to evaluate energies of DCH, CVV and WW states. 
The experimental ionization energies of the ss-DCH states have been measured for several 
molecules^iii and have been well reproduced by the CASSCF method.-ii^ Also, the CVV 
and WW energies of NH3 obtained by the CASSCF and CASCI calculations roughly coin- 
cide with the peak positions in the experimental DCH Auger spectrum.— In order to obtain 
theoretical Auger electron kinetic energy distribution, we have to calculate Auger intensi- 
ies other than the energies of the DCH, CVV and WW states. In the present work, we 
directly evaluate Auger amplitude in the Wentzel's formula, where the initial and the final 
states are represented by multi-configuratinal wave functions. As we will show in section II, 
the Auger amplitude is represented by a linear combination of two-electron integrals with 
core-hole, valence and continuum Auger electron orbitals. In the literature of molecular sin- 
gle core-hole (SCH) Auger decay, explicit evaluation of continuum Auger orbital has been 
limited to a few theoretical works,JS"— since treatment of such orbital in molecular system 
is rather difficult. Instead, several approximate methods have been introduced to avoid this 
difficulty. For instance, the molecular orbital of Auger electron is replaced by an continuum 
Auger orbital in atomic system.— The Stieljes imaging method was also employedSS, to 
represent continuum orbital in basis set. In other case, atomic population of valence 
orbitals is utilized to estimate absolute or relative Auger transition intensities,—"— in which 
continuum orbital is not evaluated at all. In this work, the intensities of DCH — )■ CVV 
and CVV — )■ WW Auger decays are approximated by using atomic population of valence 
orbitals, which was originally introduced for molecular SCH Auger decay by Mitani et al.— 
Although this approximation is not so accurate, it will be suitable to survey characteristic 
feature of Auger electron spectrum originated from molecular DCH Auger decay. 
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II. THEORETICAL METHOD 



A. Auger Intensity 

Kinetic energies of Auger electrons can be evaluated using standard quantum chemistry 
method for bound electronic state. In the present work, we employ the CASSCF and CASCI 
methods to obtain Auger electron energies. Estimation of Auger transition intensity, on the 
other hand, is not so straitforward as in evaluation of Auger electron kinetic energy. As in 
other theoretical works in molecular core-hole decay,— intensities of DCH — )■ CVV 
Auger decay and subsequent CVV — )■ WW Auger decay are evaluated by the Wentzel's 
formula,— 



Ifi = 271 



{^j\H-E\^,) =2n\t\\ (1) 



where and are the wave functions for the initial and the final state, H is the Hamil- 
tonian, E is the energy of the initial state, and t is the amplitude of the Auger transition. 
Atomic units are used in this expression. The wave function for the initial core-hole state, 
is represented by a single or multi configurational electron function with bound molec- 
ular orbitals (MOs), where is the number of electrons in the initial state. In contrast, the 
wave function for the final state is in general represented as = AJ2a^a.(pa, where $a is 

— 1 electron wave function for the final ion state a, (pa represents continuum function for 
Auger electron in channel a, and A is antisymmetlization operator. In this work, we employ 
single-channel expression of the final state in which we include only one final ion state 
$a- When we calculate an intensity of Auger transition between the initial state and the 
final ion state $a, the final state wave function is just represented as = A^a4>a under 
the single-channel expression. Further discussion on this subject can be found in Ref. 
and references therein. In the present work, we used DCH and CVV state for and 
respectively, in calculation of DCH — ?► CVV Auger intensity. Similarly, we used CVV and 
WW state for and ^a, respectively, in calculation of CVV WW Auger intensity. 

Before going into detail of actual Auger intensity calculation, we show below that con- 
venient short expressions are available for DCH CVV Auger amplitude, when the initial 
DCH state and the final CVV ion state are represented by single configurational state 
functions (CSFs) and the frozen orbital approximation is adopted. Although these short 
expressions were not used to calculate Auger intensity in the present work, we used them 



to interpret and analyze our Auger spectra. In case of SCH Auger decay, it is well known 
that the Auger amplitude t in Eq. ([1]) has convenient short expressions when the initial and 
the final ion states are represented by single CSFs.— '22i>^i^ Assuming that the doublet SCH 
state is prepared by inner-shell ionization of closed-shell molecule, the Auger amplitude for 
the singlet final ion state with two vacancies in the v and w valence MOs is, 

t = -^^[{kvlcw) + {kw\cv)] {v w) , (2) 
= {kv\cv) {v = w) . (3) 

For triplet final ion state, the amplitude is 

R 

t = Y - [{kv\cw) — {kw\cv)] . (4) 

Here c and k represent MOs for the core-hole vacancy and the continuum Auger electron, 
respectively. The expression {ij\kl) represents two-electron integral, 

{ij\kl) = / dridr2(j)i {n) (ri) —(pk (rz) (f)i (rz) , (5) 

which involves MOs (pi, (pj, (pk and (pi. In a similar manner, we can derive expression for 
amplitudes of Auger transition from DCH state to CVV state, assuming frozen orbital ap- 
proximation as well as single CSF wave functions. The results are summarized in Table 
[H We selected a specific spin-coupling scheme to represent doublet CVV states in Table 
[H where intermediate spin state, singlet or triplet, was first formed in valence electrons, 
then it was coupled with doublet core electron to form total spin state with S=l/2. In 
other word, we represent two S=l/2 wave functions as $5 = l/\/2{\vwc\ — \vwc\) and 
= 1/ {\vwc\ + \vwc\ —2\vwc\). Spin-coupling scheme is not unique when we con- 
struct total spin state from three S=l/2 particles.— For example, the wave functions $5 
and $T shown above are constructed from spin coupling scheme vw ■ c, which means that v 
and w orbitals form singlet or triplet configuration, then c orbital couples this configuration. 
Similarly, other spin coupling scheme such as v ■ wc or cv ■ w is possible. So, care must be 
taken when expressions in Table [I] are compared with other calculation. The expressions in 
Table [Hare very similar to those of the SCH Auger amplitudes, i.e., addition of two-electron 
integrals appears in the expression when valence electrons are coupled to form singlet, while 
subtraction of the integrals appears when valence electrons form triplet. From this obser- 
vation, the Auger transition from DCH state to the CVV state with singlet intermediate 



valence spin state is expected to be more intense than the transition to the CW state with 
triplet intermediate spin state, as in the case of SCH Auger decay. Although it may be 
possible to derive similar expressions for CVV — > WW Auger decay, we did not attempt 
it because there are too many combinations in orbitals. In the present work, we directly 
evaluated Eq. ([T]) using CI wave functions of CVV and WW states, as will be described 
below. Although the equations in Table [T] were not used in our actual calculation, they may 
be usefuU in other theoretical or experimental work when quick evaluation of DCH — t- CVV 
Auger intensity is required. 

When the initial and the final ion states are represented by more general multi configu- 
rational functions, the Auger amplitude is represented as 

t = ^Cy^^^{kv\cw) , (6) 

v,w 

where the summation of MOs v and w is taken over all active valence orbitals. The coef- 
ficients Cy^w depend on CI coefficients of the wave functions, determinant-CSF conversion 
coefficients, and relative ordering of MOs in determinants of the initial and final states. As 
in the case of single CSF wave functions, frozen orbital approximation was assumed in Eq. 
(E]). We employ this frozen orbital approximation because of its computational efficiency. If 
we used non-orthogonal orbital sets, considering relaxation effect, then evaluation of Auger 
amplitude would be much difficult because of long calculation time. Discussion on validity 
of the frozen orbital approximation for single core-hole Auger decay can be found in Ref. |30 
and references therein. In order to obtain Auger spectrum, we directly evaluated the Auger 
amplitudes t in Eqs. ([1]) and (|6]) in combination with multi-configurational expressions (CI 
wave functions) of the initial and the final states. Eqs.(Il]) and ([6]) can be used for any type 
of Auger transition involving one core-hole decay, two-valence hole creation and an emission 
of one Auger electron, which includes DCH— t-CW and CVV— t- WW Auger decays as well 
as SCH— i-VV Auger decay. It is important to recognize that C^,^^ in Eq. ([6]) contains the 
product of the CI coefficients of the initial and final ion state wave functions, thus only the 
CSFs with large CI coefficients contribute to the Auger amplitude t. 

The two-electron integral {kv\cw) in Eqs. ([3D,(11D and (jH]) contains MO k for the continuum 
Auger electron. We used the method described in Mitani et al.— which approximates this 
two-electron integral {kv\cw) by the population of MO w on the atom A where the core- 
hole orbital c is localized. Although this approximation is not so accurate, it can reproduce 
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qualitative feature of the experimental SCH Auger spectra of H2O and NH3,— and will be 
appropriate to understand qualitative nature of DCH Auger spectrum as well. 

B. Detail of the calculation 

The energies of DCH, CW and WW states as well as Auger intensities were evaluated 
by fixed bond length calculations. Since the lifetime of DCH and CW states are expected 
to be similar or shorter than the lifetime of SCH state, e.g., less than 10 fs for C, N or O 
elements,— the effect of the vibrational motion may be small. Experimental geometries of 
CH4, NH3 and H2CO molecules^"— were used in this work. Calculation for H2CO molecule 
was performed with C2V point group symmetry. Although CH4 and NH3 molecules belong 
to Tfi and C^^ point group symmetries, respectively, we used D2 and Cg symmetries in the 
calculations. The energies and CI coefficients of the CW and WW states were obtained 
by the CASCI calculations using the MOs prepared by the state-averaged (SA) CASSCF— 1^ 
calculations for the low-lying CW and WW states. For CH4, 60 CW and 300 WW 
states were obtained by the CASCI calculations based on the CASSCF calculations for 
the low- lying 12 CW and 10 WW states, respectively. For NH3, 100 CW states were 
calculated by the CASCI method based on the CASSCF calculation for the low-lying 15 CW 
states, and 300 WW states were obtained by the CASCI method using the MOs taken from 
the CASSCF calculation for the low-lying 11 WW states. For H2CO, 1500 Cls^^ CW, 
1500 Ols^^ CW and 3000 WW states were obtained by the CASCI calculations based 
on the CASSCF calculations for the low-lying 20 Cls"^ CW, 24 Ols-^ CW and 30 WW 
states, respectively. Note that these numbers refer to total numbers of calculated CW 
or WW states in the molecules, i.e., summation of calculated states over all irreducible 
representations. These electronic states obtained by the CASCI calculations were sufficient 
to cover whole energy range of the DCH — t- CW as well as the CW — )■ WW Auger spectra 
of CH4, NH3 and H2CO. The cc-pVTZ basis set^^ was employed for all these calculations. 
We used frozen Hartree-Fock orbitals to represent Is core orbitals, while the other orbitals 
were fully relaxed and optimized in the CASSCF calculations. As in our previous works on 
molecular DCH states,-"^ all valence electrons were distributed in the active orbital space 
composed of all available valence MOs, while occupation of core-hole orbital was explicitly 
restricted. The same set of active orbitals was used in both the CASSCF and CASCI 
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calculations. The total number of active orbitals is 8, 7 and 10 for CH4, NH3 and H2CO, 
respectively. With these active orbitals, typical number of configurations in the CVV states 
is 670, 550 and 12500 for CH4, NH3 and H2CO, respectively. For the WW states, the 
number of CSFs is 90, 100 and 3500 for CH4, NH3 and H2CO, respectively. The energies of 
the DCH states for CH4, NH3 and H2CO were also obtained by the CASSCF calculations, 
using the same basis set and active orbitals as in the calculations for the CVV and WW 
states. Our CASSCF double ionization energies (DIEs) are 650.2, 891.0, 658.0 and 1168.5 eV 
for the CH4 Cls-2, NH3 Nls-^, H2CO Cls'^, and H2CO Ols'^ ss-DCH states, respectively. 
The DIE of the H2CO Cls-^Ols-^ ts-DCH state is 846.6 eV for the singlet, and 847.0 eV for 
the triplet state. These DIEs were used to obtain the Auger kinetic energies of the DCH — )■ 
CVV Auger decays. For these CASSCF and CASCI calculations, molpro program package^ 
was mostly used. In addition, congen and scatci modules in the UK R-matrix codesAl were 
partly used to analyze the electronic states. 

In evaluation of the Auger intensities, we used the Lowdin atomic population of the MOs 
which were taken from the CASSCF calculation for the CVV states. The population ob- 
tained by the MOs for the WW states gave similar Auger intensities. Although choice 
of population, Lowdin or MuUiken, affects Auger intensities as demonstrated by Mitani et 
al.,— qualitative feature of the result does not change much. In Auger intensity evaluation 
using Eqs. ([1]) and (|6]), wave function of DCH state was approximated by a single CSF, 
since weights of the other configurations were small. For CVV and WW states, five con- 
figurations with the largest CI coefficients were taken into account in the Auger intensity 
evaluation. We performed a test calculation on CH4 DCH Auger spectra, taking into ac- 
count 30 configurations with the largest CI coefficients. The results obtained with 30 and 
5 configurations are almost the same, as shown in Fig. SI. This indicates that inclusion of 
five configurations is enough in Auger intensity evaluation. The intensities of the 1st Auger 
transitions were normalized as ^ ■ / (DCH, CWj) = 1 or 3, where the summation is 1 for 
the singlet ss- and ts-DCH states and 3 for the triplet ts-DCH state, and I (DCH, CWj) is 
the intensity for transition from the DCH state to the final state with the ith CVV ion state. 
The normahzation for the 2nd Auger intensities is Y.j I (CW^, VVVVj) = / (DCH, CW^), 
where / (CWj, WWj) is the intensity for the transition from the zth CVV state to the 
final state with the jth WW ion state. This expression for normalization of the 2nd Auger 
intensity indicates that two conditions should be satisfied at once in order to have large in- 
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tensity in the 2nd Auger decay; (1) The ith CVV state is populated enough in the 1st Auger 
decay, or, / (DCH, CVVj) is large enough, and (2) The ith CVV state and the jth WW 
ion state have enough transition moment in Eq. ([T]), or, / (CVVj, VVVVj) is large enough. 
These two conditions severely restrict possible pair of CVV and WW states which has 
noticeably large intensity in the 2nd Auger decay. 

III. RESULTS AND DISCUSSION 
A. CH4 

In Fig. [21, calculated CH4 Auger intensities are shown as a function of Auger electron 
kinetic energy, where the 1st Auger transitions from the Cls^^ DCH state to the Cls~^ 
CVV states and the 2nd Auger transitions from the Cls^^ CVV states to the WW 
states contribute the spectrum. Since the main configuration of the CH4 DCH state is 
(ls)°(2ai)^(lt2)^, the 2ai and 1^2 valence orbitals mainly participate to the Auger transi- 
tions. Main part of the DCH — t- CVV Auger spectrum extends from 270 to 290 eV, and 
the CVV — i- WW Auger spectrum extends from 210 to 240 eV. These two different Auger 
spectra are well separated in energy each other. Number of discrete Auger transitions, in- 
dicated as vertical bars in Fig. HI is 6 for the DCH CVV Auger decay and 300 for the 
CVV —7- WW Auger decay. The convoluted DCH — t- CVV Auger spectrum in Fig. [21 
obtained by convolution of the discrete Auger spectrum with Gaussian having 4.5 eV width, 
has 3 large distinct peaks, whereas the convoluted CVV — > WW Auger spectrum has 4 
peaks. The highest energy peak (~ 295 eV) in the convoluted DCH — t- CVV Auger spectrum 
corresponds to the (^2)"^ valence vacancy creation. The second highest (~ 280 eV) and the 
3rd highest (~ 270 eV) energy peaks in the DCH — t- CVV spectrum are formed by the 
(201)^^(^2)^^ and (2ai)^^ valence hole creations, respectively. The highest (~ 235 eV) and 
the lowest (~ 210 eV) energy peaks in the CVV — )■ WW Auger spectrum originate from 
the (^2)"^ and (2ai)^^ vacancy creations, as in the DCH — )■ CVV Auger case. Two peaks 
in the middle of the CVV — )■ WW Auger spectrum are formed by the (201)^^(^2)^^ va- 
cancy creation, where the initial states of the higher energy peak (~ 225 eV) have the (2ai)^ 
configurations and those of the lower energy peak (~ 220 eV) have the (2ai)^ configurations. 

The kinetic energies of the 1st and 2nd Auger electrons have correlation because these 
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two Auger decays are not independent as shown schematically in Fig. [H The upper panel of 
Fig. [3] shows 2D Auger intensity distribution as functions of the 1st and 2nd Auger electron 
kinetic energies, obtained by smoothing theoretical discrete Auger spectrum. Integration of 
this 2D Auger spectrum along the vertical axis (the 2nd Auger electron energy) gives the 
convoluted DCH — t- CVV Auger spectrum in Fig. ^ and integration along the horizontal 
direction (the 1st Auger electron energy) gives the CVV — WW Auger spectrum. In 
the upper panel of Fig. |3l several distinct and round-shaped high-intensity regions are 
recognized. The highest intensity region around the 1st Auger energy of 295 eV and the 2nd 
Auger energy of 235 eV corresponds to the (lt2)^^ valence hole creation followed by another 
(1^2)^^ valence hole creation in the successive two Auger transitions. Next to this highest 
intensity peak, there are three peaks with the 2nd highest Auger intensities. The peak 
around the 1st Auger energy of 295 eV and the 2nd Auger energy of 220 eV corresponds to 
the (1^2)"^ vacancy creation in the 1st Auger decay, followed by the (2ai)~^(lt2)~^ vacancy 
creation in the 2nd Auger decay. The peak at the 1st Auger energy of 280 eV and the 2nd 
Auger energy of 235 eV corresponds to successive formation of the (2ai)~^(lt2)~^ and the 
{lt2)~'^ valence vacancies. The peak around the 1st Auger energy of 280 eV with the 2nd 
Auger energy of 225 eV is formed by the (2ai)-^(lt2)"^ vacancy creation followed by another 
(2ai)^^(lt2)~^ vacancy creation. Note that Lablanquie et al. reported the experimental 2D 
Auger spectrum for the ss-DCH Auger decay of N2 molecule,— by measuring kinetic energies 
of the two Auger electrons. So far, such experimental 2D Auger spectrum is not available 
for the other molecules. 

Binding energies of DCH state (i^Dcn), CVV state (i^cvv) and WW state (i^vvvv) 
are related to the 1st and 2nd Auger electron kinetic energies (-EAugcri and -EAuger2) as, 
Edch = EcYv + -^Augcri and -Ecvv = -Eyvvv + -E'Auger2- Usiug these relations as well as 
experimental or theoretical value of -Edch, the 2D Auger spectrum in the upper panel of 
Fig. [3] can be converted to the 2D Auger spectrum as functions of CVV and WW binding 
energies, as shown in the lower panel of Fig. [31 In this plot, contributions of CVV and 
WW states to the Auger intensity can be easily recognized: the highest intensity peak at 
bottom left corner of the plot is related to the lowest energy CVV and WW states, for 
example. By integrating this 2D Auger spectrum along the horizontal axis (CVV energy), 
integrated one-dimensional (ID) Auger intensity as a function of WW binding energy can 
be obtained as shown in the lower panel of Fig. |H Similarly, integration along the vertical 
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axis (WW energy) gives integrated ID Auger intensity as a function of CW binding 
energy as shown in the upper panel of Fig. HI Although we need integration to extract these 
ID spectra from experimental data, theoretical ID spectra can be obtained directly from our 
calculation, without explicit integration. The original discrete spectra are shown as vertical 
bars in Fig. HJ Three distinct peaks are seen in both panels, with the highest energy peak 
corresponds to the (2ai)°(l)f:2)" configuration, the second highest peak corresponds to the 
(2ai)^(lt2)"~^ configuration and the lowest energy peak corresponds to the (2ai)^(lt2)"~^ 
configuration, where n is 6 for the CW states in the upper panel and 4 for the WW states 
in the lower panel. 

For reference, ionization energies of the calculated CW states are summarized in Table 
[TTlwith their main configurations and intensities for the 1st Auger decay. The WW states 
are not provided because there are too many states. Because of spin conservation in Eq. (Q, 
the singlet ss-DCH state of CH4 decays only to the doublet CW states. When the wave 
functions of the DCH and CW states are approximated by a single CSF, expressions of the 
ss-DCH — i- CW Auger intensities can be classified into three different types as described in 
Sec. II. A and Tabled According to this classification. Auger intensity tends to be larger for 
the CW final ion state with singlet intermediate spin state in valence electrons than that 
with triplet intermediate spin state in valence electrons. The calculated Auger intensities in 
Table [III roughly obey this trend. 

B. NH3 

In Fig. [5l calculated Auger spectrum of NII3 DCH decay is shown as a function of Auger 
electron kinetic energy, which include contributions from the 1st Auger transitions from the 
Nls~^ DCH state to the Nls^-*^ CW states and the 2nd Auger transitions from the Nls~-^ 
CW states to the WW states. For comparison, experimental NH3 DCH Auger spectrum 
of Eland et alM and theoretical Auger spectrum of NH3 SCH — )■ W (doubly valence ionized 
states) Auger decay are also shown in the figure. The SCH Auger spectrum was calculated 
by the same procedure as we used for the DCH — )• CW and CW — t- WW Auger decays. 

Main part of the theoretical DCH — )■ CW Auger spectrum extends from 380 to 420 eV, 
whereas the CW — )■ WW and SCH W Auger spectra extend from 300 to 350 eV and 
335 to 375 eV, respectively. Although the CW WW and SCH W Auger spectra 
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overlap partly, the main part of the DCH — > CVV Auger spectrum is well separated from the 
CVV —7- WW and SCH — > W spectra. Number of discrete Auger transitions, indicated 
as vertical bars in Fig. El is 14 for the DCH — )■ CVV Auger decay, 23 for the SCH — )■ VV 
Auger decay, and 995 for the CVV — ?■ WW Auger decay. Because of difference in number 
of transitions, discrete spectra of the DCH — )■ CVV and CVV — i- WW Auger decays 
look rather different in Fig. [51 However, when these spectra are convoluted by Gaussian 
function, close similarity is observed between these two Auger spectra. Also, overall shape 
of the convoluted SCH VV Auger spectrum resembles those of the convoluted DCH — )■ 
CVV and CVV — t- WW Auger spectra. Compared to the position of the convoluted SCH 
— VV Auger spectrum, the DCH — )■ CVV Auger spectrum is located about 50 eV higher in 
energy, and the CVV — ?■ WW spectrum is located about 20 eV lower in energy. Origin of 
the energy shift of the DCH — t- CVV Auger spectrum relative to the SCH — t- VV spectrum is 
mainly attributed to the Coulomb repulsion of the core holes, which increases the energy of 
the DCH state.-"^ In addition, the CVV states are stabilized by existence of two valence holes 
which weaken the Coulomb screening of the nuclear charges. This stabilization contributes to 
the shift of the DCH — t- CVV Auger spectrum into the higher energy. The shift of the CVV 
— 7> WW Auger spectrum to the lower energy can also be explained by this stabilization 
of the CVV states. Each convoluted Auger spectrum has 5 distinct peaks, which can be 
assigned in terms of the occupied valence MOs in the main configurations of the DCH and 
SCH states: the 2ai, le and 3ai orbitals. The highest energy peak has the largest intensity, 
accompanying a faint shoulder structure at higher energy side. This shoulder structure 
corresponds to the Auger transition with (3ai)^^ valence vacancy creation, whereas the 
main peak corresponds to the (le)~^(3ai)~^ vacancy creation. The second highest energy 
peak, with the second largest intensity, originates from the (le)~^ vacancy creation. The 
3rd, 4th and 5th peaks are formed by the {2ai)~^ {3ai)'^ , (2ai)~^(le)^^ and (2ai)~^ valence 
hole creations, respectively. The intensity of the 5th peak is much smaller than the other 
peaks. 

The experimental spectrum of Eland et al.,— obtained by the triple electron coincidence 
method with synchrotron radiation, is compared with our results in Fig. [5l Although 
the experimental spectrum was plotted as a function of CVV binding energy in Ref. uM, we 
converted the spectrum as a function of Auger electron kinetic energy using the experimental 
NH3 DCH binding energy of 892 eV. The large peak around 425 eV in the experimental 
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spectrum coincides well with the (3ai)~^, (le)~^(3ai)~^ and (le)~^ peaks in the calculated 
DCH —7- CVV Auger spectrum. The rise of experimental intensities around 400 eV may 
have relation with the (2ai)^^(3ai)^^ and (2ai)~^(le)~^ peaks in the calculated spectrum. 
The experimental intensity appears to increase at around 380 eV, where the calculated 
(2ai)^^ peak exists. However, association of this 380 eV peak with the calculated peak 
is not clear because of difference in intensities as well as experimental uncertainty. The 
experimental peaks around 360 - 370 eV coincide well with the calculated SCH — t- VV peaks 
with the (3ai)~^, (le)~"'^(3ai)^^ and (le)^^ vacancy creations. As noted by Eland et al.,— 
the origin of these peaks around 360 - 370 eV is secondary electrons created by the SCH 
— VV Auger electrons or photoelectrons of the SCH formation. Thus, the height of these 
peaks do not directly reflect the SCH formation cross section, which is much larger than the 
DCH formation cross section. Also, the shape of the SCH related peaks may be distorted 
from the original SCH Auger spectrum, because of the process involved in the secondary 
electron creation. We just adjusted the height of the calculated SCH Auger spectrum to 
the experimental peak around 370 eV. There is a peak around 340 eV in the experimental 
spectrum, whose location is very close to the calculated CVV — )■ WW Auger peaks with 
the (3ai)~^, (le)^"'^(3ai)^"'^ and (le)^^ vacancy creations. Although intensities are smaller, 
the (2ai)^^(3ai)^^ and (2ai)~^(le)~^ peaks in the SCH — )■ VV Auger spectrum may also 
have contribution to the experimental peak at 340 eV. Below 330 eV, association of the 
calculated CVV — i- WW Auger spectrum and the experimental spectrum is not clear. 
Although agreement with the experimental spectrum and our result looks modestly good on 
the whole, some degree of discrepancy exists, especially around the 2nd Auger component. 
Origin of the difference may be related to the approximations employed in our calculation, 
e.g., the frozen orbital approximation or neglect of the shake-up satellite DCH states. By 
considering orbital relaxation or the satellite states, the shape of the theoretical Auger 
spectrum may be changed. In addition, experimental noise from SCH signal can be another 
source of the discrepancy, which may be reduced by performing four-electron coincidence 
experiment. 

In the upper panel of Fig. [HI 2D Auger spectrum is shown as functions of the 1st and the 
2nd Auger kinetic energies. Compared to the 2D Auger intensity distribution for the CH4 
DCH decay in Fig. [3l high-intensity regions look less distinct in the NH3 Auger spectrum. 
For both the 1st and 2nd Auger electrons, the higher kinetic energy corresponds to vacancy 
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creation in the 3ai and le orbitals, and the lower kinetic energy corresponds to vacancy 
creation in the 2ai orbital. Thus, the highest-intensity region surrounded by the 1st Auger 
energies of 410 - 420 eV and the 2nd Auger energies of 340 - 350 eV mainly corresponds to 
the (le3ai)~^ valence hole creation followed by another (le3ai)^^ valence hole creation in 
the sequential two-step Auger decays. The second highest intensity regions are located next 
to the highest intensity region, one surrounded by the 1st Auger energies of 410 - 420 eV and 
the 2nd Auger energies of 320 - 335 eV, and the other surrounded by the 1st Auger energies 
of 390 - 405 eV and the 2nd Auger energies of 340 - 350 eV. The former corresponds to the 
(le3ai)~^ valence hole creation followed by the (2ai)^^(le3ai)^^ valence hole creation, and 
the later corresponds to the (2ai)~^(le3ai)^^ valence hole creation followed by the (le3ai)~^ 
valence hole creation. 

As in the CH4 case, the 2D Auger spectrum in the upper panel of Fig. [6] can be converted 
to the 2D intensity as functions of CVV and WW binding energies as shown in the lower 
panel of Fig. El This plot shows clearly that the highest intensity region corresponds to 
the low-lying CVV and WW states. Integration of the 2D Auger spectrum in the lower 
panel of Fig. |6] along the vertical or the horizontal axis gives the integrated ID Auger 
spectra as shown in Fig. [71 The integrated Auger intensity as a function of CVV binding 
energy in the upper panel of Fig. [71 has 3 large peaks which correspond to the CVV states 
with (2ai)^, (2ai)^ and (2ai)° configurations. The integrated Auger intensity as a function 
of WW binding energy in the lower panel has more complicated structure compared to 
the intensity in the upper panel. We can still make assignment of these peaks in terms 
of electronic configuration, e.g., the lowest energy peak corresponds to the configuration 
(2ai)^(le)^(3ai)°, the second lowest peak to the configuration (2ai)^(le)^(3ai)^, while the 
peaks in the middle of the spectrum are related to the (2ai)^ type configuration. Compared 
to the integrated ID spectra of CH4 WW states, the integrated ID spectrum of NH3 
WW states is distributed broader and looks less structured. This difference is attributed 
to the different number of occupied valence orbitals in CH4 (two orbitals) and NH3 (three 
orbitals). This suggests that it may be hard for large molecule to recognize peak structure 
in integrated ID spectra. 

For reference, ionization energies of the calculated CVV states are listed in Table [TTTl with 
their main configurations and the intensities for the DCH — )■ CVV Auger decay. As in the 
CH4 case, the DCH — )• CVV Auger intensity tends to be larger for the CVV final ion states 
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with singlet intermediate valence spin than the ion states with triplet intermediate valence 
spin. 

C. H2CO 

In Fig. [HI the Auger electron kinetic energy distribution of H2CO is shown for the 
sequential Auger decays of the DCH and CVV states. In contrast to NII3 and CII4, H2CO 
has four different DCH states: two ss-DCH states ( singlet Cls~^ and Ols^^ ) and two ts- 
DCH states ( singlet and triplet Cls^^Ols^^ ). We consider the case where these four DCH 
states are created in XFEL experiment. In order to simplify the situation, we just assumed 
that the formation probabilities are determined by statistical ratio: 1:3 for singlet and triplet 
DCH states. This approximation has been well known and discussed in the literate of inner- 
shell photoionization.— 1^ The formation probabilities of the singlet ts-DCH and the ss-DCH 
states are assumed to be the same, because they are expected to be formed through sequential 
two-photon two-electron ionization. We also calculated the Auger spectra expected in SR 
experiment, where only the ss-DCH states contribute the Auger spectra because of very low 
formation probabilities of the ts-DCH states. The result is shown in Fig. S3. 

The lower panel of Fig. [8] shows the Auger spectra associated with the Cls core-hole 
decays, which include the Auger transition from the Cls~^ ss-DCH state to the Cls^^ CVV 
states, the transitions from the Cls^^Ols"^ ts-DCH states to the Ols^^ CVV states, and 
the transitions from the Cls^^ CVV states to the WW states. Similarly, the upper panel 
of Fig. [8] shows the Auger spectra associated with the 01s core-hole decays, including the 
Auger transition from the 01s~^ ss-DCH state to the Ols^^ CVV states, the transitions 
from the Cls^^Ols^^ ts-DCH states to the Cls^^ CVV states, and the transitions from the 
Ols"^ CVV states to the WW states. Since the 01s ^ CVV states are generated by the 
Auger decays of the Ols^^ ss-DCH and the Cls-^Ols^i ts-DCH states, the Ols^^ CVV ^ 
WW Auger spectrum has two different origins. Likewise, the Auger decays of the Cls~^ 
ss-DCH and the Cls-^Ols-^ ts-DCH states contribute to the Cls^^ CVV ^ WW Auger 
spectrum. The DCH — t- CVV and CVV — ?■ WW Auger decays of H2CO involve mainly 
6 valence orbitals, i.e., three ai, one 61 and two 62 MOs, because the main valence electron 
configuration of the DCH states is represented as (3ai)^(4ai)^(162)^(5ai)^(16i)^(262)^- Since 
this number of orbitals is larger than those in the CH4 and NH3 Auger decays, the discrete 
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H2CO Auger intensities in Fig. [8] are more densely distributed in comparison with the 
intensities for the CH4 and NH3 Auger decays. Also, peaks in the convoluted H2CO Auger 
spectra are not so distinct as in the spectra of the CH4 and NH3 Auger decays. 

As shown in the lower panel of Fig. [8], the main part of the Cls~^ ss-DCH — t- Cls^^ 
CVV Auger spectrum extends from 270 to 310 eV, with weak intensities around 250 eV. 
The Cls^^Ols"^ ts-DCH — )■ Ols^^ CVV Auger spectrum, including singlet and triplet con- 
tributions, extends from 210 to 250 eV, with weak intensities around 190 eV. The spectrum 
of the Cls^^ CVV —7- WW Auger decays extends from 180 to 260 eV. As can be seen in 
the figure, the Auger spectrum of the Cls~^ ss-DCH Cls^^ CVV decay is well separated 
from the other spectra. However, the main parts of the Cls^^Ols^^ ts-DCH — )■ Ols^^ CVV 
and the Cls~^ CVV — ?■ WW Auger spectra overlap each other. The main part of the 
Ols"^ ss-DCH — )■ Ols^^ CVV Auger spectrum in the upper panel of Fig. [8] extends from 
510 to 570 eV, whereas the Cls~^01s"^ ts-DCH — )• Cls~^ CVV Auger spectrum extends 
from 460 to 500 eV. As in the case of the Cls^^ core-hole decays, the spectrum of the 01s~^ 
CVV -> WW Auger decay overlaps with those of the Cls^^Ols-i ts-DCH ^ Cls^^ CVV 
Auger spectrum. Our results suggest that the Auger decays of the ss-DCH states may be 
easily identified in experiment. However, distinction between the ts-DCH CVV Auger 
spectra and the CVV — t- WW Auger spectra may be difficult. 

The convoluted Cls~^ ss-DCH — )• Cls^^ CVV Auger spectrum has two large peaks at 300 
eV and 275 eV, and one faint peak at 250 eV. The peak at 300 eV is formed by two-electron 
vacancy creation in the 5ai, I61, I62 and 262 valence MOs, (5ail6il622&2)^^, while the peak 
at 275 eV corresponds to the (3ai)~^(5ail6i 162262)"^ valence hole creation. The structure 
between the 300 eV peak and the 275 eV peak involves vacancy creation in the 4ai MO as 
well. The faint peak at 250 eV is formed by the (3ai)~^ vacancy creation. Note that these 
assignments are not as obvious as in the cases of the NH3 and CH4 Auger spectra. Similarly, 
we can relate the peaks in the other convoluted Auger spectra to the (5ail6i 162262)"^; 
(3ai)^^(5ail6il62262)^^ and (3ai)~^ vacancy creations, with contribution of the 4ai hole 
creation to the structure between the highest energy and the second highest energy peaks. 
For example, the peak around 560 eV in the 01s~^ ss-DCH — )■ Ols^^ CVV Auger spectrum 
and the peaks around 240 eV and 480-490 eV in the Cls-^Ols-^ ts-DCH CVV Auger 
spectra correspond to the (5ail6il62262)^^ valence hole creation. Also, the peaks around 
530-540 eV in the 01s~^ ss-DCH — )• Ols^^ CVV Auger spectrum and the peaks around 210- 
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220 eV and 460-470 eV in the Cls"^01s-^ ts-DCH CVV Auger spectra are formed from 
the (3ai)^^(5ail6il62262)^^ vacancy creation. The origin of the peaks in the convoluted 
CVV — )■ WW Auger spectra were not inspected because these peaks contain too many- 
discrete transitions. Since the shapes of the convoluted CVV — ?■ WW spectra are roughly 
similar to the DCH — t- CVV Auger spectra, origin of the peaks may be the same as in the 
DCH CVV Auger spectra. 

In the upper panel of Fig. [9l 2D Auger intensity distribution for the Cls~^ ss-DCH 
Cls^^ CVV Auger decay and the subsequent Cls^^ CVV — t- WW Auger decay is shown 
as functions of the 1st and the 2nd Auger electron kinetic energies. Similarly, 2D Auger 
spectrum for the Ols"^ ss-DCH ^ Ols-^ CVV and the Ols-^ CVV ^ WW Auger decays 
is shown in the lower panel of Fig. |9l These 2D Auger spectra for the ss-DCH decays may 
be available in experiment using 4-electron coincidence technique, since the Auger spectra 
originated from the ss-DCH — t- CVV Auger decays are separated from the other Auger 
spectra as seen in Fig. [HI Compared to the peaks in the 2D spectra of the CH4 and NH3 
DCH Auger decays, the peaks in the H2CO ss-DCH decays look less clear. The highest 
intensity regions are located around the 1st Auger energies of 295 - 305 eV and the 2nd 
Auger energies of 220 - 240 eV for the Cls~^ case, and around the 1st Auger energies of 
550 - 565 eV and the 2nd Auger energies of 470 - 490 eV for the 01s~^ case. Location of 
these highest intensity regions indicates that the (5ai 161162262)^^ valence hole creation and 
subsequent (5ail6il622&2)^^ valence hole creation are the most probable in the successive 
two Auger decays. The 2D Auger spectra in Fig. [9] are converted to the 2D intensities 
as functions of CVV and WW binding energies as shown in Fig. [TOl Based on these 
spectra, the highest-intensity peaks can be related to the low-lying CVV and WW states. 
In addition, we can recognize that some CVV states with moderately high energies, 385 eV 
in the Cls^^ case and 640 eV in the 01s~^ case, contribute to the higher-intensity peaks 
in Fig. [ini These 2D intensities as functions of CVV and WW binding energies may be 
available in experiment, if the binding energies of the ss-DCH states as well as the 2D Auger 
intensities as functions of the 1st and 2nd Auger electron energies are measured. 

In the two-step Auger decays of the Cls"^01s~^ ts-DCH states, two different Auger 
decay pathways exist: (1) the 1st Auger transitions from the ts-DCH states to the 01s~^ 
CVV states, and subsequent Auger transitions from the Ols^^ CVV states to the WW 
states, (2) the 1st Auger transitions from the ts-DCH states to the Cls^^ CVV states, and 
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subsequent Auger transitions from the Cls~^ CW states to the VWV states. In the case 

(1) , the 1st Auger electron has an energy in the region of the Cls core-hole decay and 
the 2nd Auger energy is in the region of the 01s core-hole decay. In contrast, in the case 

(2) , the 1st Auger electron has an energy in the region of the 01s core-hole decay and 
the 2nd Auger energy is in the region of the Cls core-hole decay. The 2D Auger intensity 
distributions originated from these two different Auger decay pathways are shown in the 
panel (a) and (b) of Fig. [11] as functions of the 1st and 2nd Auger electron energies. It 
may be difficult in experiment to determine the order of two successive Auger decays, or 
in other word, to make distinction between these two Auger decay pathways. Thus, two 
Auger intensities of the case (1) and (2) will not be observed separately, but summation of 
these two intensities will be measured as shown in the panel (c) of Fig. [TTJ The highest 
intensity region is located around the Cls Auger energy of 230 - 240 eV and the 01s Auger 
energy of 490 eV in the panel (c). As in the cases of the ss-DCH — )■ CVV Auger decays, 
this region roughly corresponds to 4-electron valence vacancy creation in the 5ai, Ibi, I62 
and 262 orbitals in the sequential two Auger decays. Although other peaks and structures 
appear in the 2D intensity distribution, it is difficult to make clear assignment because two 
different Auger intensities are added in the panel (c). In contrast to the 2D Auger spectra of 
the ss-DCH decays, the 2D Auger spectrum of the ts-DCH decays in the panel (c) cannot be 
converted to the intensity as functions of CVV and WW binding energies, because such 
conversion requires individual 2D Auger spectra as in the panels (a) and (b). The integrated 
ID Auger spectra as a function of VWV binding energy, as in the lower panels of Figs. H] 
and [71 can be obtained from the 2D spectra in the panels (a) and (b), using the relation 
-E'vvvv = -f'DCH — -f'Augeri " -£'Auger2- Becausc this relation is preserved by exchange of -EAugeri 
and -EAuger2, the ID Auger spectrum can be directly obtained from the 2D spectrum in the 
panel (c). Thus, if experimental value of -Edch is available, information on WW binding 
energies may be obtained in experimental measurement on the ts-DCH Auger decays. 

For reference, ionization energies of the low-lying Cls^^ CVV states are shown in Table 
IIVI with their main configurations and Auger intensities. Also, the low-lying 01s~^ CVV 
states are listed in Table [Vl As in the ss-DCH — > CVV Auger decays of NH3 and CH4, 
the singlet DCH states decay only to the doublet CVV states. In addition to these doublet 
states, the triplet Cls^^Ols^^ ts-DCH states decay to the quartet CVV states as well. In 
the DCH —7- CVV Auger decays of NH3 and CH4, the Auger intensity tends to be larger for 
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the CVV final ion state with singlet intermediate valence spin than the state with triplet 
intermediate valence spin. However, this tendency is not obvious in H2CO. 

D. Discussion 

Our calculation on NH3 DCH Auger spectrum has roughly reproduced the experimental 
Auger spectrum of Eland et al.,^- and thus similar agreements are expected for CH4 and 
H2CO DCH Auger spectra. The convoluted Auger spectra contain several large distinct 
peaks, which can be related to two hole creations in valence MOs. In case of the integrated 
ID Auger spectra, we can relate their peaks to the electronic configurations of the CVV 
or WW states. Although this kind of assignment can be efficiently performed for small 
molecules such as CH4 and NH3, it may be difficult for larger molecule with many occupied 
valence MOs. 

Since there are many CVV and WW states in CH4,NH3 and H2CO, possible number of 
transitions in the CVV — ?■ WW Auger decays might be huge, and thus distributions of the 
Auger spectra might be much broader and less structured compared to those of the DCH 
— i- CVV Auger decays. In contrast to this intuitive expectation, the distributions of the 
calculated CVV — ?■ WW Auger spectra have similar extent as in the DCH — > CVV Auger 
spectra, as can be seen in Figs. [2], [5] and [H In addition, we can see several distinct peaks in 
the convoluted CVV — ?■ WW Auger spectra, which is especially pronounced in the CH4 
case. Our results can be understood in the following manner. When we evaluated the DCH 
— 7- CVV Auger intensities using Eqs. (1) and (6), we approximated the DCH state by a single 
CSF function. This means that the DCH — > CVV Auger intensity is noticeably large only 
when the CI wave function of the CVV state is dominated by a configuration obtained by 
removing two valence electrons from and adding one core electron to the DCH configuration, 
because only such a configuration has non-zero Auger amplitude in Eq. (1). Similarly, the 
CVV — 7- WW Auger intensity is noticeably large only when (1) the CVV state has large 
intensity in the DCH — )■ CVV Auger decay and (2) the CI wave function of the WW 
state is dominated by a configuration obtained by removing two valence electrons from and 
adding one core electron to the main configuraton of the CVV state. This also means that 
the final WW ion states are all dominated by electronic configurations constructed from 
the occupied valence MOs of the DCH state. The above conditions (1) and (2) severely 
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restrict possible pairs of CVV and WW states which have large contribution to the CW 
— > WW Auger decay. The distribution of the CVV — )■ WW Auger spectra cannot 
be as broad as intuitively expected from the distribution of the CVV and WW binding 
energies, because of the selection conditions stated above. These conditions also indicate 
that the Auger electron kinetic energy in the CVV — i- WW Auger decay is related to the 
two-electron valence hole creation in the occupied MOs of the DCH state. In other word, 
the CVV — > WW Auger intensity can be noticeably large only around the specific energy 
region related to the two-electron valence hole creation, although this property may be less 
pronounced as number of occupied MOs increases. 

The total intensity of the H2CO ts-DCH Cls^^ CVV Auger transitions, obtained by 
integrating the spectrum by the Auger kinetic energy, is larger than that of the H2CO ts- 
DCH — > Ols^^ CVV Auger transitions, as seen in Fig. El This difference indicates that 
the occupied valence MOs of the ts-DCH states have larger overlap with the oxygen site 
than that with the carbon site. If individual, separated experimental spectra are available 
for ts-DCH ^ CVV, ss-DCH ^ CVV and CVV -> WW Auger transitions of H2CO, and 
if assignment of peaks is available in terms of the occupied valence MOs, we can compare 
intensities of peaks having the same MO assignment in different Auger transitions. Such 
comparison may be useful in understanding relaxation of the valence MOs and overlap of 
MOs with different atomic sites. Similar analysis for ss- and ts-DCH Auger decays of other 
molecules will be interesting, if experimental measurement is possible. 

In this work, the Wentzel's ansatz was employed to estimate intensities of the 1st Auger 
transition from DCH to CVV states, and of the 2nd Auger transition from CVV to WW 
states. This means that we assume a "three-step" model in which DCH state creation by 
photoionization, the 1st Auger decay, and the 2nd Auger decay are all treated independently, 
as in the two-step model of single core-hole Auger decay.— 1^""— This assumption is valid when 
the lifetimes of the intermediate states, DCH and CVV states in this case, are sufficiently 
long that the interaction between photoelectron and the 1st Auger electron as well as the 
interaction between the 1st Auger electron and the 2nd Auger electron are negligibly small. 
According to the multiconfiguration Dirac-Fock calculation by Chen on atomic DCH state,— 
the Auger transition rates of DCH states for the low-Z elements are much larger than the 
Auger rates of SCH states; in the case of neon, the Auger rate per K hole of the DCH state 
is about 42 % larger than the rate of the SCH state. The Auger lifetimes of the CH4, NH3 
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and H2CO DCH states may also be short, but such values are not known so far. Thus, it 
will be interesting as well as necessary to verify the validity of the "three-step" assumption 
in future, for example, by performing calculation or experimental measurement on lifetime 
of DCH and CVV states, or by measuring angular distribution of Auger electrons in the 
molecular-frame as in the SCH Auger process.— 

Following Mitani et al.,^^ we approximated two-electron integrals in Eq. (|6]), which 
involve core-hole, valence and Auger continuum orbitals, by the Lowdin population of the 
valence MO in the integral. When an Auger transition amplitude is dominated by a pair 
of CSFs, one in the initial state and the other in the final state, this approximation may 
be valid as demonstrated by Mitani et al.— However, when more than one pair of CSFs 
equally contributes to an Auger amplitude, this approximation can be less accurate because 
it does not necessarily preserve proper relative phases of two-electron integrals in the Auger 
amplitude. In our calculation. Auger amplitudes were dominated by one pair of CSFs in 
most cases and thus this approximation is expected to be valid. Still, refinement of this 
approximation is desirable. 

In this paper, we have studied basic properties of molecular DCH Auger decay by es- 
timating Auger spectra originated from the normal Auger decays of the molecular DCH 
states. Angular distribution of the Auger electron was not studied in this work, though the 
experimental measurement has already been performed for the Auger decay of the N2 ss- 
DCH state.- Angular distribution of the Auger electron may contain additional information 
on valence holes created in DCH or CVV Auger decays, and will be an interesting subject 
for theoretical study in future. In addition to the normal Auger decay of molecular DCH 
state, the Auger decays of doubly core excited resonance state and shake-up satellite state 
have been studied experiment ally^iii and may deserve further theoretical study as well. 

At this point of time, experimental Auger spectra are only available for the Auger decay 
of ss-DCH state.->^»^ However, the XFEL experiment on CO molecule has been performed 
at LCLS and the Auger spectra of the ts-DCH decay are being analyzed.— Thus, we can 
expect further experimental information in this subject in the near future. 
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IV. SUMMARY 



We have performed theoretical investigation on Auger decay of molecular DCH state, for 
which experimental information becomes available for several molecules recently. Assuming 
sequential two-step Auger transitions from DCH state to CVV states and from CVV states to 
WW states, Auger electron kinetic energy distributions were estimated for CH4, NH3 and 
H2CO molecules based on the CASSCF and CASCI calculations. The Auger spectra of the 
CH4 and NH3 ss-DCH Auger decays contain two well separated components: one from the 
1st Auger transition from the DCH state to the CVV states and the other from the 2nd Auger 
transition from the CVV states to the WW final ion states. Our result roughly agrees with 
the experimental Auger spectra of NH3 ss-DCH decay,— but experimental spectrum with 
better energy resolution is desired for precise comparison. The calculated Auger spectrum 
of H2CO DCH decay has more complicated structure compared to the spectra of CH4 and 
NH3 DCH decays, due to existence of the ts-DCH Auger decay. In the Auger spectrum of 
H2CO DCH decay, the components originated from the ss-DCH CVV Auger decays are 
well separated from the rest of the spectrum. However, the components originated from the 
ts-DCH — !■ CVV Auger decays, and the components from the CVV WW Auger decays 
overlap each other, making separation of the spectra difficult. The 2D Auger spectrum may 
be helpful in resolving this difficulty. We hope our calculation and analysis on the results, 
especially the H2CO DCH Auger spectrum, will be useful in interpreting experimental result 
expected in the near future. 
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1 An example of two-step Auger decays originated from the ss-DCH state of 
NH3. Thin black arrows represent displacements of electrons during the Auger 
decays. In the 1st Auger decay, the CVV state with one vacancy in the 3ai 
orbital and another vacancy in the le orbital is produced. Then in the 2nd 
Auger transition, this CVV state decays to produce the WW state with 
additional vacancies in the le and 2ai orbitals 

2 Calculated intensities of the Auger decays of the CH4 Cls^^ DCH state 
(DCH^CW) and the Cls-^ CVV states (CVV^VVVV), as a function of 
Auger electron kinetic energy. The vertical full lines represent the discrete 
Auger spectrum obtained by the CASCI wave functions with the frozen orbital 
approximation. The dashed lines are obtained by convoluting the discrete 
Auger intensities with Gaussian function having 4.5 eV width 

3 (Upper panel) Two-dimensional (2D) intensity of the two-step Auger decays 
of the CH4 Cls~^ ss-DCH state, as functions of kinetic energies of the 1st 
and the 2nd Auger electrons. The 1st electrons are emitted by the Auger 
transition from the Cls~^ DCH state to the Cls~^ CVV states, and the 
2nd electrons are emitted by the transition from the Cls^^ CVV states to 
the WW states. This 2D spectrum represents a probability of finding two 
Auger electrons at specific pair of energies. (Lower panel) 2D Auger intensity 
as functions of CVV and WW binding energies, converted from the 2D 
spectrum in the upper panel. These intensities were obtained by smoothing 
calculated discrete Auger intensities, as in the convoluted intensities in Fig. 
[2J Unit of the intensity is arbitrary. 

4 Integrated ID Auger intensities as a function of CVV binding energy (upper 
panel) and WW binding energy (lower panel). These spectra can be ob- 
tained by integrating the 2D Auger spectrum in the lower panel of Fig. [3] by 
WW energy (upper panel) or CVV energy (lower panel). The vertical bars 
represent original theoretical intensities which can be directly obtained from 
our calculation. The other details are the same as in Fig. [2j 
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5 Calculated and experimental Auger intensities of the NH3 core-hole decays as 
a function of Auger electron kinetic energy. DCH— t-CVV: the transition from 
the Nls-2 DCH state to the Nls-^ CVV states, CVV^VVVV: the transition 
from the Nls^^ CVV states to the WW states, SCH^VV: the Auger decay 
of the Nls"^ SCH state. Expt.:experimental Auger spectrum of Eland et al.— . 
The other details are the same as in Fig. [2l 

6 (Upper panel) 2D intensity of the two-step Auger decays of the NH3 Nls~^ 
DCH state, as functions of kinetic energies of the 1st and the 2nd Auger 
electrons. (Lower panel) 2D Auger intensity as functions of CVV and WW 
binding energies, converted from the 2D spectrum in the upper panel. The 
other details are the same as in Fig. [31 

7 Integrated ID Auger intensities as a function of CVV binding energy (upper 
panel) and WW binding energy (lower panel), obtained by integrating the 
2D Auger spectrum in the lower panel of Fig. O The other details are the 
same as in Fig. HI 

8 (Upper panel) The Auger intensities for the 01s core- hole decays of H2CO 
molecule, which include the Auger decays of the H2CO Ols^^ ss-DCH state, 
the Cls^^Ols^^ ts-DCH states, and the Ols^^ CVV states. (Lower panel) 
The Auger intensities for the Cls core-hole decays of H2CO molecule, which 
include the Auger decays of the H2CO Cls^^ ss-DCH state, the Cls^^Ols"^ 
ts-DCH states, and the Cls~^ CVV states. The other details are the same as 
in Fig. H 

9 (Upper panel) 2D intensity of the two-step Auger decays of the H2CO Cls~^ 
ss-DCH state, as functions of kinetic energies of the 1st and the 2nd Auger 
electrons. (Lower panel) 2D intensity of the two-step Auger decays of the 
H2CO 01s~^ ss-DCH state, as functions of kinetic energies of the 1st and the 
2nd Auger electrons. The other details are the same as in Fig. [31 
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10 (Upper panel) 2D Auger intensity of the H2CO Cls~^ ss-DCH decay as func- 
tions of CVV and WW binding energies, converted from the 2D spectrum 
in the upper panel of Fig. [HI (Lower panel) 2D Auger intensity of the H2CO 
Ols^^ ss-DCH decay as functions of CVV and WW binding energies, con- 
verted from the 2D spectrum in the lower panel of Fig. O The other details 
are the same as in Fig. [31 

11 2D Auger spectra of the two-step Auger decays of the H2CO Cls~^01s~^ ts- 
DCH states. Panel (a): 2D spectrum with the 1st Auger transitions from the 
Cls^^Ols^^ state to the Cls^^ CVV states and the 2nd transitions from the 
Cls-^ CVV states to the WW states. Panel (b): 2D spectrum with the 1st 
Auger transitions from the Cls~^01s^^ state to the 01s~^ CVV states and 
the 2nd transitions from the Ols^^ CVV states to the WW states. Panel 
(c): 2D Auger intensity as functions of Cls and 01s Auger electron energies, 
obtained by adding the intensities of two different Auger decay pathways in 
the panel (a) and (b) 

51 Effect of number of CI coefficients included in the calculation. We performed 
the same calculation as in Fig. 2, but with 30 configurations with the largest 
CI coefficients. The result is very similar to the CH4 DCH Auger spectra 
in Fig. 2 which were calculated with 5 configurations with the largest CI 
coefficients. The details of the figure is the same as in the Fig. 2 

52 Effect of Gaussian width used for convolution. The convoluted CII4 DCH 
Auger spectra with different FWHM width are shown. Our result suggests 
that FWHM has strong influence on the shape of the Auger spectra, which 
may be important when we compare calculated spectrum with experimental 
result. The details of the flgure is the same as in Fig. 2 

53 H2CO DCH ^ CVV and CVV ^ WW Auger spectra expected in SR 
experiment. In Fig. 8, we show H2CO DCH Auger spectra expected for XFEL 
experiment. In case of SR experiment, formation efficiencies of the ts-DCH 
states are expected to be low, 1/100 of the ss-DCH states. Thus, we can ignore 
the ts-DCH contributions in simulating Auger spectra in SR experiment. In 
this figure, the Auger spectra without the ts-DCH contribution is shown. The 
details of the figure is the same as in Fig. 8 
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FIG. 1. An example of two-step Auger decays originated from the ss-DCH state of NH3. Thin 
black arrows represent displacements of electrons during the Auger decays. In the 1st Auger 
decay, the CVV state with one vacancy in the 3oi orbital and another vacancy in the le orbital is 
produced. Then in the 2nd Auger transition, this CW state decays to produce the YWY state 
with additional vacancies in the le and 2ai orbitals. 
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FIG. 2. Calculated intensities of the Auger decays of the CH4 Cls"^ DCH state (DCH^CVV) 
and the Cls^"*^ CVV states (CVV— t-VVVV), as a function of Auger electron kinetic energy. The 
vertical full lines represent the discrete Auger spectrum obtained by the CASCI wave functions 
with the frozen orbital approximation. The dashed lines are obtained by convoluting the discrete 
Auger intensities with Gaussian function having 4.5 eV width. 
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FIG. 3. (Upper panel) Two-dimensional (2D) intensity of the two-step Auger decays of the CH4 
Cls^^ ss-DCH state, as functions of kinetic energies of the 1st and the 2nd Auger electrons. The 1st 
electrons are emitted by the Auger transition from the Cls~^ DCH state to the Cls~^ CVV states, 
and the 2nd electrons are emitted by the transition from the Cls~^ CVV states to the WW 
states. This 2D spectrum represents a probability of finding two Auger electrons at specific pair 
of energies. (Lower panel) 2D Auger intensity as functions of CVV and WW binding energies, 
converted from the 2D spectrum in the upper panel. These intensities were obtained by smoothing 
calculated discrete Auger intensities, as in the convoluted intensities in Fig. [2j Unit of the intensity 
is arbitrary. 
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FIG. 4. Integrated ID Auger intensities as a function of CVV binding energy (upper panel) and 
WW binding energy (lower panel). These spectra can be obtained by integrating the 2D Auger 
spectrum in the lower panel of Fig. [3] by WW energy (upper panel) or CVV energy (lower panel). 
The vertical bars represent original theoretical intensities which can be directly obtained from our 
calculation. The other details are the same as in Fig. [21 
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FIG. 5. Calculated and experimental Auger intensities of the NH3 core-hole decays as a function 
of Auger electron kinetic energy. DCH^CVV: the transition from the Nls"^ DCH state to the 
Nls"^ CVV states, CVV^VVVV: the transition from the Nls"^ CVV states to the WW states, 
SCH— T-VV: the Auger decay of the Nls^^ SCH state. Expt.: experimental Auger spectrum of Eland 
et al.— . The other details are the same as in Fig. [2l 
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FIG. 6. (Upper panel) 2D intensity of the two-step Auger decays of the NH3 Nls~^ DCH state, 
as functions of kinetic energies of the 1st and the 2nd Auger electrons. (Lower panel) 2D Auger 
intensity as functions of CVV and WW binding energies, converted from the 2D spectrum in the 
upper panel. The other details are the same as in Fig. O 
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FIG. 7. Integrated ID Auger intensities as a function of CVV binding energy (upper panel) and 
WW binding energy (lower panel), obtained by integrating the 2D Auger spectrum in the lower 
panel of Fig. [6l The other details are the same as in Fig. SI 
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FIG. 8. (Upper panel) The Auger intensities for the 01s core-hole decays of H2CO molecule, which 
include the Auger decays of the H2CO Ols^^ ss-DCH state, the Cls^^Ols^^ ts-DCH states, and 
the 01s~^ CVV states. (Lower panel) The Auger intensities for the Cls core-hole decays of H2CO 
molecule, which include the Auger decays of the H2CO Cls~^ ss-DCH state, the Cls^^Ols^^ 
ts-DCH states, and the Cls^^ CVV states. The other details are the same as in Fig. [21 
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FIG. 9. (Upper panel) 2D intensity of the two-step Auger decays of the H2CO Cls~^ ss-DCH 
state, as functions of kinetic energies of the 1st and the 2nd Auger electrons. (Lower panel) 2D 
intensity of the two-step Auger decays of the H2CO 01s~^ ss-DCH state, as functions of kinetic 
energies of the 1st and the 2nd Auger electrons. The other details are the same as in Fig. [3j 
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FIG. 10. (Upper panel) 2D Auger intensity of the H2CO Cls~^ ss-DCH decay as functions of CVV 
and WW binding energies, converted from the 2D spectrum in the upper panel of Fig. [9l (Lower 
panel) 2D Auger intensity of the H2CO 01s~^ ss-DCH decay as functions of CVV and WW 
binding energies, converted from the 2D spectrum in the lower panel of Fig. [9j The other details 
are the same as in Fig. [3j 
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FIG. 11. 2D Auger spectra of the two-step Auger decays of the H2CO Cls^ Ols^ ts-DCH states. 
Panel (a): 2D spectrum with the 1st Auger transitions from the Cls^-'^Ols^"'^ state to the Cls~^ 
CVV states and the 2nd transitions from the Cls"^ CVV states to the WW states. Panel (b): 
2D spectrum with the 1st Auger transitions from the Cls~"'^01s~^ state to the Ols^"*^ CVV states 
and the 2nd transitions from the Ols^^ CVV states to the WW states. Panel (c): 2D Auger 
intensity as functions of Cls and 01s Auger electron energies, obtained by adding the intensities 
of two different Auger decay pathways in the panel (a) and (b). 
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FIG. SI. Effect of number of CI coefficients included in the calculation. We performed the same 
calculation as in Fig. 2, but with 30 configurations with the largest CI coefficients. The result is 
very similar to the CH4 DCH Auger spectra in Fig. 2 which were calculated with 5 configurations 
with the largest CI coefficients. The details of the figure is the same as in the Fig. 2. 
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FIG. S2. Effect of Gaussian width used for convolution. The convoluted CH4 DCH Auger spectra 
with different FWHM width are shown. Our result suggests that FWHM has strong influence on 
the shape of the Auger spectra, which may be important when we compare calculated spectrum 
with experimental result. The details of the figure is the same as in Fig. 2. 
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FIG. S3. H2CO DCH ^ CVV and CVV ^ WW Auger spectra expected in SR experiment. 
In Fig. 8, we show II2CO DCH Auger spectra expected for XFEL experiment. In case of SR 
experiment, formation efficiencies of the ts-DCH states are expected to be fow, 1/100 of the ss- 
DCH states. Thus, we can ignore the ts-DCH contributions in simulating Auger spectra in SR 
experiment. In this figure, the Auger spectra without the ts-DCH contribution is shown. The 
details of the figure is the same as in Fig. 8. 
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TABLE I. Expressions of Auger decay amplitude t for transitions from ss-DCH, singlet ts-DCH and triplet ts-DCH states to CW state, 
with assumptions of single configurational wave functions and frozen orbitals. "CVV S=l/2 (S)" represents doublet CVV state with singlet 
intermediate spin coupling of valence electrons, and "CVV S=l/2 (T)" represents doublet CVV state with triplet intermediate spin coupling 
of valence electrons, v and w refer to orbitals involved in the valence hole creation, c stands for inner-shell orbital involved in the core-hole 
decay, and k represents continuum orbital of Auger electron. 





ss-DCH 


ts-DCH (singlet) 


ts-DCH (triplet) 




{kv\cw) + (kwlcv) {v / w) 


[{kv\cw) + {kw\cv)] {v / w) 


[{kv\cw) + {kw\cv)] {v / w) 


CW S=l/2 (S) 


\pl {kv\c'w) (y = w) 


(kvlcw) {v = w) 


{kv\cw) {v = w) 


CVV S=l/2 (T) 


V3[{kv\cw) — {kw\cv)] 


•y/I [{kv\cw) — {kw\cv)\ 


[{kv\cw) — {kw\cv)] 


CVV S=3/2 






[{kv\cw) — {kw\cv)\ 



TABLE II. Representative CW states of CH4 molecule. Energy refers to the ionization energy 
with respect to the neutral ground state of CH4. Intermediate spin means spin coupling in valence 
electrons, where S and T represent singlet and triplet, respectively. 



State 


Energy (eV) 


Main configuration (Intermediate spin) 


Intensity (arb.) 




352.662 


(ls)i(lt2)^ (T) 


- 




354.362 


{is)Hit2r (s) 


0.000 




354.809 


(l.)i(li2)^ (T) 


0.000 




355.358 


M\lt2)' (S) 


0.338 


'A, 


356.676 


{ls)\lt2)' (S) 


0.165 




361.996 


(l.)i(2oi)Hlt2)^ (T) 






365.469 


(ls)H2a,y{lt2r (T) 


0.075 


^T2 


370.181 


(ls)i(2ai)HU2f (S) 


0.372 


'A, 


378. 9G1 


(l.s)i(2«i)" (S) 


0.018 
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TABLE III. Representative CVV states of NH3 molecule. The other details are the same as in 
Table ini 



State 


Energy (eV) 


Main configuration (Intermediate spin) 


Intensity (arb.) 




470.431 


(l.)i(le)3(3ai)i (T) 


- 


Ml 


470.810 


(l.)i(3ai)0 (S) 


0.078 




474.026 


(l.)Hle)3(3ai)i (S) 


0.168 




475.207 


(l.)Hle)3(3ai)i (T) 


0.154 




476.080 


(ls)i(le)2 (T) 


- 




479.240 


(l.)i(le)2 (T) 


0.000 




479.320 


(l.)i(le)2 (S) 


0.134 


'A, 


480.593 


(ls)i(le)2 (S) 


0.140 


^A, 


482.519 


(ls)i(2ai)i(3ai)i (T) 


- 


% 


487.973 


(l.)i(2ai)i(le)i (T) 




Ml 


488.317 


(ls)i(2ai)i(3ai)i (T) 


0.049 


Ml 


491.578 


(ls)i(2ai)i(3ai)i (S) 


0.066 


2e 


492.381 


(l.)i(2ai)i(le)i (T) 


0.056 




497.398 


{ls)H2a,y{le)' (S) 


0.128 


'A, 


508.638 


(l.)i(2ai)0 (S) 


0.020 
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TABLE IV. Low-lying Cls"^ CVV states of H2CO molecule. Cls-^ Cls-^Ols-^ (S) and Cls-^Ols"^ (T) are the initial states of the 1st 
Auger transitions, and represent the Cls~^ ss-DCH state, the singlet and triplet Cls~^01s~^ ts-DCH states, respectively. The other details 
are the same as in Table HTl 



State 


Energy (eV) 


Main configuration (Litermediate spin) 


Cls-2 


Intensity (arb.) 
Cls-^Ols-i (S) 


Cls-iQls-i (T) 




347.931 


(2ai)i(262)° (S) 


0.007 


0.029 


0.083 


^A, 


351.299 


{2a,)Hlb,)H2b2y (T) 






0.000 




351.646 


{2ai)Hlh)H2h)' (T) 


0.000 


0.002 


0.010 


'A, 


352.141 


{2a,)Hlb2)H2b2y (T) 






0.024 




352.408 


(2ai)i(262)H5ai)^ (T) 






0.001 




352.537 


(2ai)i(262)M5ai)i (T) 


0.000 


0.001 


0.000 


'A, 


353.532 


(201)1(162)^262)^ (T) 


0.002 


0.042 


0.009 


'A2 


353.091 


(2gi)1(16i)H2&2)^ (S) 


0.063 


0.068 


0.200 


'B, 


354.603 


(2ai)i(262)H5ai)i (S) 


0.053 


0.084 


0.198 


'A, 


355.407 


{2ai)Hlb2)H2b2y (S) 


0.097 


0.005 


0.063 



TABLE V. Low-lying Ols"^ CVV states of H2CO molecule. The other details are the same as in Table IIVI 



otate 


Energy (eV) 


Main configuration (Intermediate spin) 


Uls 


Intensity (arb.) 
Uis Uis (bj 


Uis Uis ( i j 


2 A 


595.700 


(laij (202)'' (Sj 


U.UOl 


0.040 


0.080 


4 A 


5yo. / 11 


(laij (l&ij (2ft2j (Ij 






n n /I o 
0.049 




599.230 


(101)^162)^262)^ (T) 






0.037 




599.395 


(lai)i(16i)n2?>2)^ (S) 


0.000 


0.030 


0.024 




599.876 


(lai)i(262)i(5ai)i (T) 






0.006 


^i?2 


600.951 


(lai)i(262)H5ai)i (T) 


0.004 


0.001 


0.009 




60L788 


(lai)i(162)n262)^ (T) 


0.016 


0.013 


0.001 




602.113 


(lai)i(262)H5ai)i (S) 


0.010 


0.004 


0.037 


2^2 


602.145 


(lai)i(16i)i(262)i (T) 


0.077 


0.001 


0.034 


'B^ 


602.876 


(lai)i(262)i(4ai)i (T) 






0.025 




603.166 


(lai)i(162)H2?>2)^ (S) 


0.041 


0.007 


0.061 



